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ABSTRACT 


We present the results of high resolution (R>30,000) optical and near-IR spec¬ 
troscopic monitoring observations of HBC722, a recent FU Orionis object that 
underwent an accretion burst in 2010. We observed HBC722 in optical/near-IR 
with the BOES, HET-HRS, and IGRINS spectrographs, at various points in the 
outburst. We found atomic lines with strongly blueshifted absorption features 
or P Cygni prohles, both evidence of a wind driven by the accretion. Some lines 
show a broad double-peaked absorption feature, evidence of disk rotation. How¬ 
ever, the wind-driven and disk-driven spectroscopic features are anti-correlated 
in time; the disk features became strong as the wind features disappeared. This 
anti-correlation might indicate that the rebuilding of the inner disk was inter¬ 
rupted by the wind pressure during the hrst two years. The Half-Width at 
Half-Depth (HWHD) of the double-peaked prohles decreases with wavelength, 
indicative of the Keplerian rotation; the optical spectra with the disk feature are 
htted by a G5 template stellar spectrum convolved with a rotation velocity of 
70 km s“^ while the near-IR disk features are htted by a K5 template stellar 
spectrum convolved with a rotation velocity of 50 km s“^. Therefore, the op¬ 
tical and near-IR spectra seem to trace the disk at 39 and 76 Rq, respectively. 
We ht a power-law temperature distribution in the disk, hnding an index of 0.8, 
comparable to optically thick accretion disk models. 

Subject headings: Stars; formation — Stars; FU Orionis — Optical; spectroscopy — 
Individual; HBC 722 
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1. Introduction 


During star formation, accretion onto protostars through disks seems to occur 
episodically (Audard et ah 2014 and references therein), and FU Orionis-type objects 
(hereafter, FUors) have been proposed as prominent examples of burst accreting protostars. 


HBC 722, a new me mber of the F 


latter half of 2010 fjSemkov et ah 


Jor gr oup, brightened by about 5 magnitudes (in V) in the 


20101 1. HBC 722 reac hed a relative ma ximum brightness 


in September 2010 and then faded until May 2011 (e.g. 


Sung et ah 


2013 1. Late in 2011 the 


brightness began to increase gradually and reached its second maximum in the hrst half of 
2013, after which it remained at approximately constant level. The object shows sho rt-term 


brightness variations with several day periods flGreen et ah 


2013a : 


Baek et al 


201511 . The 


broad observational coverage and clear delineation of several epochs within the optical light 
curve make HBC 722 a good case study of the accretion process and its effect on low mass 
star formation. 


Followup observations were not restricted to optical variability. 


Miller et al. 


(1201111 


characterized the source with TripleSpec in th e near-IR, and obs erved a typical spectrum 


for a T Tauri star, with a high accretion rate. 


Green et al 


Herschel spectroscopy (at submillimeter wavelengths) and 


(1201 111 observe c 


Dunham et al. 


the source with 


( 201211 observed 


the source with the SMA at millimeter wavelengths; both concluded that the source was 


a relatively evolved young star, lacking a signihcant ex tended envelope, a nc 


accret ion rate (few x 10 ® AfQyr ^) by FUor standards (iKospal et al.ll201111 . 


a modest 


Liebhart et al. 


(1201411 observed HBG 722 with Ghandra, and found that the X-ray increased during the 
burst, concluding that dust-free gas had been ejected by the star sometime after the early 


^Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project 
of the University of Texas at Austin, the Pennsylvania State University, Stanford University, 
Ludwig-Maximilians-Universitat Miinchen, and Georg-August-Universitat Gottingen. 
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outburst. Taken together, this suggests that the source was a low mass Class II (T Tauri) 
star pre-outburst, and it underwent a sudden accretion event that initiated a jet sometime 
around or after the initial optical burst. 


The triggering event for FUors is not fully understood; in particular it is unclear 
whether the burst proceeds from the very innermost regions of the circumstellar disk, or 
whether it begins as a perturbation a t larger disk radii following the formation and release 


of a disk instability flZhu et al 


20071) . Optical/near-IR variability can trace inner disk 


activity but is difficult to dis 


flickering; e.g.. 


Kenvon et al 


ingu ish from the activity on or near the stellar surface (e.g. 


20001) . It is only with high resolution spectroscopy, covering 


the earliest stages of burst, that we can distinguish variation in the disk from variations 
near the star and outflow/jet. 

In this paper, we report monitoring spectroscopic observations of HBC 722 with 
high resolution spectrometers in optical and near-IR. We use higher cadence but lower 
signal-to-noise spectroscopy to guide our less frequent high signal-to-noise spectroscopic 
observations. We hnd a highly blueshifted wind component in absorption in the profiles 
of various lines inclnding Ha and Ca II A8498, which show a P Cygni prohle. We also 


observe a broad doub 


e-peaked absorption fe ature, which is an evidence of the disk rotation 


flHartmann fc Kenvon 


1996 


Zhu et al. 


20091 ). in the neutral metal lines both at the optical 


and near-IR wavelengths. These double-peaked features imply a rotating Keplerian inner 
disk (within < 1 AU) for which we can constrain the current temperature prohle. 
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2. Observations 

2.1. HET - HRS 


We observed HBC 72 2 with the High Re solution Spectrograph flTull 


Hobby-Eberly Telescope flRamsev et ah 


1998h of the 


1998h at McDonlad Observatory on the nights of 


07 December 2010, 31 May 2011, 03 August 2012 and 30 August 2012, 30 May 2013 and 
08 August 2013. These spectra were taken using the “600g5822” HRS conhguration. This 
provided spectral coverage from 479 to 681 nm at a resolving power of i? = A/AA = 30, 000. 
The target was imaged onto a 2 arcsec hber, and two additional 2 arcsec hbers were used 
to gather light from the sky immediately adjacent to the target. The CCD detectors were 
bined 2x2 pixels to give four binned pixels per resolution element. On each night, four 
separate spectra, each of 900 second duration, were obtained, for a total exposure time of 
60 minutes. Data reductions were performed with the IRAF echelle package. The four 
separate spectral images of HBC 722 were added together, using a ‘cosmic-ray’ rejection 
algorithm in order to remove noise spikes due to radiation hits on the CCDs. Spectra 
of HBC 722 and each of the two sky hbers were subtracted separately. The spectra from 
the two sky hbers were averaged, and that spectrum was subtracted from the HBC 722 
spectrum. The scaling of the sky spectrum was determined empirically to fully remove the 
Earth night-sky emission lines. 


2.2. BOAO - BOES 


The Bohyunsan Optical Echelle Spectrograph (BOES) is a high resolution echelle 
spectrograph attached to the 1.8-m optical telescope at Bohyuns an Optica l Astr onomy 


Observatory (BOAO). HBC 722 has been observed using BOES ( Lee et al.lboiD since 


November 2010 (see Table 1) with resolving power {R = ^) of 30,000 using the 300 
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/im fiber. The observed wavelength range was 3600A - lOSOOA, covering the fnll optical 
spectrnm. Two other FUors, FU Orionis and V1057 Cyg, have been observed with BOES 
on 26 Jannary 2011 and 11 September 2012, respectively. The spectra of standard stars 
used in Figure 1 were also observed with BOES. 

The observational data were reduced by the IRAF (Image Reduction and Analysis 
Facility) echelle package. Each aperture from the spectral images was extracted using a 
master flatfield image. Using the flathelding process, we corrected the interference fringes 
and pixel-to-pixel variations of the spectrum images. A ThAr lamp spectrum was used for 
wavelength calibration. 


2.3. HJST - IGRINS 


HBC722 was observed with the Immersion Grating Infrared Spectrograph (IGRINS) 
installed at the 2.7 m Harlan J. Smith Telescope (HJST) at the McDonald ob servatory o n 


20 November, 2014. IGRINS covers the full H and K bands with R=40.000 (Park et ah 

n 

20141 ). The spectra have been reduced by the IGRINS pipeline q A standard AO star (HIP 


103108) was observed immediately after the observation of HBG722 for the telluric line 
correction. 


3. Results 

In Figure [H we compare the lines with a strong blue-shifted absorption feature with 
wings greater than ~150 km s“^, showing the prohle variation throughout the burst. 
Spectra observed from burst (2010) until 03 August 2012 show this feature, after which it 


^https://github.com/igrins/plp 
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disappears over the course of a single month. 

We observe a P Cygni (wind/outflow) prohle in Ha and Ca II A8498. In Ha, we also 
observe a redshifted emission wing as well as a blueshifted absorption at velocities as large 
as ~400 km s“^. The Ca H profile is only covered by BOES and does not appear as wide, 
but we cannot rule out a similar profile due to lower S/N in the BOES spectra. In Figure [H 
the HET-HRS spectra have much higher signal-to-noise ratios compared to the BOES 
spectra so that we present only one BOES line, Ca H at 8498 A, which is not covered by 
the HET-HRS. 

If we examine the Na I D lines observed in May and August 2013 (Figure [T] 4th row, 
panels 5 and 6), the line is composed of a narrow deep component, typical of late F or G 
type stellar spectra, and a broad shallow component, which we attribute to disk rotation. 
We compare this with a Gl-type stellar spectrum (HD 188650; plotted in red) that matches 
the narrow deep component seen in the August 2013 spectra of HBG722. 

On first inspection, we also see a highly blueshifted shallow broad component in H/5 
and Fe H 5018. However, these features do not seem related to the blueshifted wind seen in 
Ha. These additional blueshifted and redshifted absorption components appear in H/3 and 
Fe H A5018 (rows 1 and 2 of Figured]). Additional redshifted components only are noted in 
Mg I A5184 (row 3 of Figure [I])- 

These components show a broad flattened structure that is likely associated with other 
lines broadened by the rotation of the Keplerian disk. These features are well matched with 
a G5 stellar template spectrum convolved with a rotational velocity (vsini) of 70 km s“^ 
as presented in blue together with the HBG722 spectra observed in May 2013. This disk 
rotation is also evident in various neutral metal lines (Ga I, Fe I, Li I) as seen in Figured] 
In these neutral lines, the high velocity blueshifted absorption seen in Ha is not present. 


It has long been suggested that a double-peaked broad a 


jsorption feature in an FUo r, 


like that seen in Figure [21 is clear evidence of disk rotation flHartmann &: Kenvon 


Zhu et al. 


1996; 


20091) . Numerous metal lines in our spectra show thi s double-peaked prohle; we 


list them along with measured HWHD (jPetrov fc Herbie 


2008 ) in Table 2. The average 


of HWHDs for optical lines is about 79 km s“^, and their standard deviation is 7.6 km 
s“^. It is not clear whether the scatter between individual lines inferred HWHDs may be 
signihcant. To calculate HWHD, we must determine the continuum level, so the uncertainty 
of the continuum level as well as the different S/N ratios for individual lines can cause the 
observed scatter. 


We fit these neutral metal line profiles by convolving a template G5 stellar spectrum 
with the rotational profile below. 


4>{Av) = 


Av 


- 1/2 


( 1 ) 


where Vmax = vsini . Most of the neutral metal lines are well-fitted by a G5-type stellar 
spectrum convolved with a rotational velocity of 70 km s“^ (shown in blue in the boxes for 
the spectra of HBG722 observed in May 2013, Figure [2]). 


Many of these prohles gradually acquire a central narrow absorption feature after 
May 2013 (Figure |2]); we attribute this feature to the central object, which has begun to 
contribute as the stellar brightness increased. The stellar lines of Ga I A6122, Fe I A6142, 
and Li I A6708 are very deep in August 2013. The linewidth of the central feature is 
reasonably similar to the G1 template, but the disk features are similar to the G5 template 
convolved with a rotational profile. One notable exception is Li A6707 which is not observed 
in the G5 template but appears in HBG 722. The reverse is true for Ga I A6439, which 
appears in the G5 template but not in HBG 722. 


We compare the spectra in early and late burst to FU Orionis and V1057 Gyg, 
two classical FUors, in Figure [3l Early in outburst, HBG 722 looks very similar to 
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FU Ori currently, with a strong blueshifted wind-driven feature and dou ble-peaked 


absorption profiles. FU Ori lacks a substanti al circumstellar envelope (e.g. 


2012; 

Green et al. 

201 3 J: 

Audard et al. 

2014) 


Dunham et ah 


20141) as does HBC722. However, the accretion 


luminosity of FU Ori has remained high (at 10-100 times higher accretion rate than 
HBC 722) for 75 years, whereas HBC 722 appears to be peaking after only a few years. 
V1057 Cyg (at ~ 40 years post-outburst) is surrounded by considerably more nebulosity 
than FU Ori or HBC 722, but still shows deep blueshifted absorption features. 


It has been argued for V1057 Cyg fjPetrov et al 


19981) and FU Orionis flPetrov fc Herbig 


20081) that the lack of difference in the width of the optical lines indicates that the 


boxy/double-peaked profiles originate from the central star only. In contrast, we attribute 
only the narrow central absorption to the central object, and the “wings” of the features 
to the disk. The detection of a separate stellar absorption in HBC 722 is possible because 
HBC 722 is a much less luminous system (5-12 Lq) than FU Ori or V1057 Cyg (100 - 250 
Lq) with a correspondingly lower inferred accretion rate. If these features are disk-driven 
in the more luminous objects, a similar narrow absorption may appear as they fade. 


Because we observed HBC 722 over multiple epochs, this dataset is the first opportunity 
to see the line profiles change as the burst evolves. The double-peaked feature in most lines 
was weak before August 2012, when the wind was still strong and dominated the absorption 
prohles. The broad blueshifted (wind) component almost disappeared in August of 2012; 
at the same time, the double-peaked profiles emerged and stabilized. 


We also observed a single epoch of HBC 722 with the high-resolution IGRINS 
spectrograph in 2014 (Figure H]). We compare spectra of HBC 722 (black) with a K5-type 
stellar spectrum (HD 44537) convolved with a rotational profile of vsini = 50 km s“^(blue). 
The narrow stellar component is clearly seen at the line center of the CO overtone 


transitions. 
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Figure [5] presents the HWHD of absorption spectra as a function of wavelength. In 
order to increase the signal-to-noise ratio, we averaged the last three epochs of optical 
spectra (30 August 2012, 30 May 2013, and 08 August 2013). The IGRINS spectra were 
observed only once. The HWHD decreases with wavelength, particularly in the IR. This is 
consistent with Keplerian disk rotation; longer wavelengths trace the cooler outer part of 
the disk where the rotation velocity is smaller. The Pearson correlation coefficient between 
HWHD and wavelength is -0.69 (p-value=0.0008), which is within 95 % confidence interval. 


4. Discussion 


Stars form by the accretion of material from disks to protostars, and this accretion 


process often triggers a disk wind. We o 


or as a P Cygni profile in optical lines flHartmann fc Kenvon 


Dserve the wind as a 


jluesh ifted absorption feature 


19961) . Our high resolution 


optical spectra clearly show a wind feature in HBC 722, which appeared as soon as the 
object was observed after its optical burst (2010). The wind velocity peaked at about 400 
km s“^, and sustained for ~2 years. It disappeared almost completely during summer 2012, 
although the brightness continued to rise into the first half of 2013 and stayed at a constant 
level after then. At the highest observed velocity, the wind could have traveled up to ~170 
AU during those 2 years. 


The second set of features are the blueshifted and redshifted boxy/double-peaked 
absorption profiles, which we attribute to a self-luminous disk rotating at Keplerian speeds. 
This interpretation is supported by the decrease in HWHD of the double-peaked profile 
from optical to near-IR (Figure |5]), and by the presence of a central narrow absorption 
feature that can be attributed to the central object. The optical lines have both greater 
rotational velocities (~ 70 km s“^ compared with 50 km s“^ in the IR) and are matched 
by an earlier spectral type template (G5) than the IR lines (K5), as presented with blue 
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spectra in Figure [2] and Figure 01 


If HBC 722 has a Keplerian disk and its central protostellar mass is 1 Mq, the 
double-peaked optical spectral trace the disk rotating at 39 Rq while the near-IR spectra 
traces the disk rotating at 76 Rq. This indicates that the temperature of the disk decreases 
from ~5,200 K at the radius of 39 Rq to ~3000 K at the radius of 76 Rq. If the power-law 
temperature distribution is adopted, this decrease of temperature corresponds to the 
power-law index of 0.8, which is very similar to | for the optically thick disk with a steady 


accretion (iHartmann fc Kenvon 


19961 ). Both temperatures are greatly enhanced over a 


quiescent T Tauri-type disk at corresponding radii, indicating that the disk has remained 
hot despite the wind turnoff in 2012. 


The bolometric luminosity of the system reached and stayed at its second maximum 
for at least 2 years beyond the wind turnoff time, which was itself 2 years after the initial 
rise (which lasted only 3 months). This suggests the initial optical burst was driven by a 
phenomenon with a few month timescale (perhaps a free fall time or a few orbital times 
at the disk inner edge), while the driving source of the wind was a longer timescale event, 
perhaps an accretion stream from further out in the disk. 


One interesting phenomenon is that the disk feature started to appear when the wind 
feature disappeared. This phenomenon might hint at the a rebuilding timescale for the 
inner disk. After a burst accretion event, the inner disk can be rebuilt by the material 
infalling from the outer disk or an extended envelope, depending on the evolutionary stage 
of the outbursting object. However, wind pressure could prevent the material from rehlling 
immediately, maintaining the inner boundary of the disk at a relatively large radius. Once 
the wind turns off, the disk material can move inward once again to rebuild the inner 
hot disk, which is then visible at optical wavelengths once more. High spatial resolution 
observations of the inner disk of HBC 722 from telescopes such as ALMA might provide 





strong constraints on disk accretion and rebuilding processes. 


5. Conclusions 

We have observed the burst of HBC 722 in unprecedented detail with high resolution 
optical and near-IR spectroscopy over multiple epochs. The optical line prohle of Ho 
reveals a blueshifted 400 km s“^ wind immediately after burst. As the system luminosity 
(and therefore accretion rate) of the system increased from 2010 - 2012, the wind remained 
strong. In 2012, the wind diminished and double-peaked rotating disk prohles (70 km s“^, 
spectral type G5) appeared in the neutral metal lines. In 2013, a narrow central absorption, 
likely the central star, appeared in the neutral metal lines. In 2014, the IR line prohles 
(including CO overtone features) were observed to be double-peaked but with a lower 
rotational velocity (50 km s“^) and later spectral type (K5), suggesting that they arise 
farther out in the disk than the optical lines. Combined with high resolution interferometry 
and additional optical/IR spectral coverage as the burst evolves, this dataset will reveal new 
constraints on the accretion and outhow process in the inner disks of low mass protostars. 
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Foundation (NRF) funded by the Ministry of Education of Korea. J.-E. L. is very grateful to 
the department of Astronomy, University of Texas at Austin for the hospitality provided to 
her from August 2013 to July 2014. JDC thanks Aditi Raye Allen for early data reduction 
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Table 1. Observing log 


Telescope 

Instrument 

(wavelength coverage) 

resolving power (R) 

target 

Date 

Exposure time 

[sec] 

HET 

HRS 

30,000 

HBC 722 

2010 Dec 07 

3600 


(4790 -6810 A) 


HBC 722 

2011 May 31 

3600 




HBC 722 

2012 Aug 08 

3600 




HBC 722 

2012 Aug 30 

3600 




HBC 722 

2013 Aug 08 

3600 

BOAO 

BOES 

30,000 

HBC 722 

2010 Nov 26 

3600 


(3600 -10500 A) 


HBC 722 

2010 Dec 11 

3600 




HBC 722 

2010 Dec 23 

3600 




HBC 722 

2010 Dec 29 

3600 




HBC 722 

2011 Jan 20 

3600 




HBC 722 

2011 Apr 20 

3600 




HBC 722 

2011 May 18 

3600 




HBC 722 

2011 Jun 16 

3600 




HBC 722 

2011 Sep 06 

3600 




HBC 722 

2011 Sep 25 

3600 




HBC 722 

2012 May 21 

3600 




HBC 722 

2012 Sep 10 

3600 




HBC 722 

2013 Apr 21 

3600 




EU Ori 

2011 Jan 26 

3600 




V1057 Cyg 

2012 Sep 11 

3600 




HD 52497 

2011 Sep 25 

300 




HD 188650 

2011 Sep 22 

600 

HJST 

IGRINS 

40,000 

HBC 722 

2014 Nov 20 

2400 


(1.47-1.81 /im, 1.95-2.48 fim) 


HIP 103108 

2014 Nov 20 

1200 




HD 44537 

2014 Nov 17 

8 








Table 2. 


List of Double-Peaked Lines 


A [A] 

Ion 

HWHD [km s-i] 

4957.60 

Fe I 

77.28 

5098.70 

Fe I 

98.41 

5302.30 

Fe I 

76.84 

5455.61 

Fe I 

82.09 

5594.47 

Ca I 

93.47 

5857.45 

Ca I 

72.66 

6056.01 

Fe I 

70.65 

6102.72 

Ca I 

80.83 

6122.22 

Ca I 

82.70 

6141.73 

Fe I 

73.03 

6191.56 

Fe I 

82.64 

6411.65 

Fe I 

75.94 

6439.08 

Ca I 

78.91 

6449.81 

Ca I 

83.24 

6462.57 

Ca I 

69.86 

6677.99 

Fe I 

71.87 

6707.89 

Li I 

72.98 

6717.68 

Ca I 

83.12 

16755 

A1 I 

63.59 

21789 

Ti I 

55.60 
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HRS 



-400 0 400 -400 0 400 -400 0 400 -400 0 400 -400 0 400 -400 0 400 

Velocity [km/s] 


BOES 

26NOV2010 25Sep2011 21May2012 10Sep2012 2lApr2013 



Fig. 1.— The time variation of optical lines including the blueshifted (wind-tracing) absorption 
feature. HET-HRS observations are presented in the upper panel, with five different lines shown 
as five rows, and six different observation epochs as six columns. The BOES observations are 
presented in the lower panel in five observation epochs, which are close to the dates for HET-HRS 
observations. The red spectra in the boxes for August 2013 is the spectra of the G1 type star, HD 
188650, which matches the best the deep narrow component of the spectra. The blue spectra in the 
boxes for May 2013 (April 2013 for BOES) is the spectra of the G5 type star, HD 52497, convolved 
with the rotational profile of vsini = 70 km s“^. The shallow broad component of each spectrum 
and the broadened adjacent features can be htted well with these convolved spectra. 
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Fig. 2.— Time variation of the double-peaked absorption features in various lines, evidence 
of disk rotation. Only the HET-HRS spectra are presented because the BOES spectra have 
much lower signal-to-noise ratio. The dates in the columns are the same as those in Figured] 
The blue and red spectra are also the same as those in Figured! 
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Fig. 3.— Comparison of HBC722 in early outburst vs. FU Ori (top row), and later in 
outburst after the apparent wind-turnoff vs. V1057 Cyg (bottom row). The HBC722 ob¬ 
servations are from HET-HRS, with the exception of the Ca II line which was observed by 
BOES. All observations of FU Ori and V1057 Cyg are from BOES. 
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Fig. 4.— IGRINS spectra (presented in black) of neutral metal lines and CO overtone 
transitions that show a double-peaked prohle. The blue spectrum is from a template K5- 
type (HD 44537) stellar spectrum convolved with a rotational velocity of 50 km s“^. The 
original stellar spectrum of HD 44537 is shown in red. 
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Fig. 5.— The HWHD of spectra with a double-peaked absorption feature (see Table 2) 
as a function of wavelength. We observe a decrease in HWHD with wavelength, consis¬ 
tent with Keplerian disk rotation. The Pearson correlation coefficient between HWHD and 
wavelength is -0.69 (p-value=0.0008), and the Kendall’s tau correlation coefficient is -0.32 
(p-value=0.05). The relation presented at the upper right is the result of linear regression 
(the solid line). 




